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We argue that a decrease of the chiral scalar meson mass is responsible for re-creation of 'I'' from 
J/^' in ultrarelativistic nucleus-nucleus collisions. This causes the charmonium yields to freeze out 
at temperatures close to the chiral symmetry restoration temperature Tc- As a result ^''/^' may 
serve as a thermometer for Tc itself. Results in a detailed reaction model support the conjecture. 
They show good agreement with recent data of NA38 and NA50 for J/'^ and ^' production in S 
on U and Pb on Pb collisions. 



Recent measurements by the NA50 collaboration at 
CERN show a so-called anomalous suppression of J/^ 
yields in non-peripheral Pb(158AGeV) on Pb collisions 
0. The results have been interpreted as a hint for 
Quark-gluon plasma (QGP) formation in these reactions 
[||- §]. The QGP can "ionize" a J/"^ analogously to 
the photo-effect j|] or even screen it out of existence 0. 
Other explanations, dissociation in coherent gluon fields 
[ 1 and collisions with secondary hadrons (comovers) 
[ 11| have also been put forward. Generally speaking, the 
discussion is centered about the question whether the 
observed suppression of charmonia is an initial or a final- 
state effect. The issue cannot be considered as settled 
yet. The arguments expressed e.g. in |2|- Q that the 
different trends in S+U versus Pb+Pb suggest the ap- 
pearence of a strong energy density-dependent destruc- 
tion mechanism for charmonium states could be tested by 
repeating the experiment at lower beam energies. Here 
we assume that the charmonium states are indeed de- 
stroyed very early |p^ . 

Naively, the fate of the different charmonium states is 
independent from each other and is governed by their 
size and binding energy. The location of the 2s excita- 
tion 5*' (3686) is barely below the threshold to the DD 
continuum while the J/^* with a mass of 3097 McV is 
deeply bound. It is natural to assume that the ^E"' is 
more easily destroyed. Indeed, NA38 has found that the 

/'^ ratio goes down continuously in S on U collisions 
with increasing centrality. However, this ratio is leveling- 
off at around 4% in non-peripheral Pb on Pb reactions 
according to the preliminary NA50 measurements. Ap- 
proximately the same value is seen in the most central S 
on U collisions. This value is roughly a factor of 4 lower 
than in pp and pA reactions. 

Does the constancy of the 4''/^' ratio result from an 
accidental fine-tuning of the parameters which govern the 
survival rates of each species? Here we answer otherwise, 
vp'/vji is leveling off, because observed ^E*' are re-created 
from J/^. If so, 5"'/^ will therefore no longer decrease 
in more violent collisions (at RHIC and LHC). The value 
of the /'^ ratio reflects on the temperature T of the 
medium provided the equilibration mechanism is suffi- 
ciently robust. In equilibrium, ^''/vE' changes with T ex- 
ponentially as exp (-(M' - M)/T) {M'/Mf^. In fact. 



/'^ is very sensitive to small changes of T, because the 
mass difference M' — M is much larger than the relevant 
temperatures. The measured /'^ ratio then suggests a 
freeze-out temperature around 170 MeV. Such a temper- 
ature is close to the expected transition temperature Tc 
for chiral symmetry restoration. 

What may the nature of the interactions which equi- 
librate the \E''/\[' ratio be? In the vacuum decays 
predominantly into a J /'^ and two pions in s-wave. The 
decay width is tiny (141 KeV), and therefore the inverse 
process in the pion gas is simply too slow. Based on 
the matrix element for the decay process we can cal- 
culate charmonium (de-) excitation probabilities in col- 
lisions with pions. However, even those are not large 
enough on the time scale of nuclear collisions. When 
^I^' J/^'TTTT decay was studied first in the seventies 
Schwinger et al. advocated an interpretation of these de- 
cays through the intermediary of a scalar resonance . 
The transition matrix element is suppressed in the vac- 
uum due to the magnitude of the scalar resonance mass. 
The presumably large transition rate between J/vt and 
in nuclear collisions may then be a consequence of 
a dropping scalar mass. The emergence of a low mass 
sigma meson which becomes degenerate with the pions 
at Tc is responsible for rapid transitions between 5"' and 
J/^P around this temperature. The transition rates be- 
come very small at slightly lower temperatures due to the 
sensitivity of the transition rates to the sigma mass. This 
means in turn that the \E''/^' ratio in heavy ion collisions 
may serve as a rather clean measure of Tc itself. 

In the remainder of this Letter we estimate the tran- 
sition rates between J/vP and in thermal matter close 
to Tc- Subsequently, we will assess the main corrections 
to the equilibrium limit. Finally, we implement the T de- 
pendent transition rates into a simple model for nucleus- 
nucleus collisions and compare the calculated J/vE* and 
^ff' yields to the prehminary NA38 and NA50 data. 

Let us consider first charmonium transitions in a heat 
bath based on the linear sigma model. At Tc the quark 
condensate vanishes and chiral symmetry gets restored. 
Seen from the hadronic side, the sigma meson mass de- 
creases with temperature T, because the sigma state has 
to become degenerate with the pion. The pion mass is 
expected not to change much with T. Its smallness is 



1 



protected by the spontaneous break-down of chiral sym- 
metry. It may be more suprising that a sigma model 
describes the relevant physics if one approaches Tc from 
above |^5|, QCD matter at very large tempera- 

tures is supposedly composed of a weakly coupled gas 
of colored quarks and gluons. However, at temperatures 
Tc < T < 2Tc, strong nonperturbative phenomena in the 
QGP may be present, since the QCD couphng constant g 
is of order one. The scalar-pseudoscalar modes are rather 
special as far as chiral symmetry is concerned. Conden- 
sation is about to occur in the a channel while the pions 
will emerge as the Goldstone bosons of the broken phase 
at slightly lower temperatures. Calculations based on the 
instanton model support the idea that correlations due 
to strong attractive forces and perhaps bound states are 
formed in the QGP at T ~ |3- 
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FIG. 1. a{'i/n 'P'tt) as a function of the out-going pion 
energy with the vp' at rest (left side). The straight (dashed) 
line represents the calculation in the symmetric phase (vac- 
uum with rUs = rrig). The right side shows the rate constant 
for ^' transitions into a J/^ due to collisions and decays 
(straight line with dropping s mass, dashed line with m^). 
The dotted line represents just the contribution from 'I'' de- 
cay in the chiral scenario. 

Brown and Cahn pointed out that the two most impor- 
tant observations - isotropic tt emission in their common 
rest frame and strong enhancement at large invariant tttt 
masses - are reproduced by the linear sigma model at 
tree level . In this model a and tt form a four- vector 
a ~ (cr, 7?) which enters into the interaction terms of the 
Lagrange density 



(1) 



A suitably chosen potential U in eq. (0) gives rise to the 
spontaneous break-down of the symmetry in which the 
isoscalar field acquires a nonvanishing expectation value, 
in the vacuum (cr)o=/w. 92 MeV is the pion decay 
constant. The effective potential U cannot be reliably 
calculated from QCD. The standard 'Mexican hat' poten- 
tial was introduced by Gell-Mann and Levy mainly 
for its simplicity and renormalization property. We are 
going to utilize the 'freedom' in the choice of U to study 



how a variation of Tc affects J/^* and ^' in AA collisions. 
In the broken phase the effective potential at finite T is 
minimized by some nonzero expectation value / of the a 
field. With the field redefinitions a = (/ + s) cos6' and 
— {f + s)tt sin 6 (with tt a vector of unit length) all bur- 
den of the self-interactions induced by U is laid on the 
scalar-isoscalar field s. s represents the fluctuations of 
the order parameter. The following interaction vertices 
which are determined from the Lagrange density enter 
into our tree-level calculation: 



V*'*00 = 2ce'^e'' 



VsTTTT = (2piP2 + ml) , Vss-^yr = -J^PlP2- 

The Cartesian isospin indices of pions have been sup- 
pressed. The coupling constant gM' of decay equals 
2c/. The field represents the s meson field in the bro- 
ken phase and an arbirary a component in the symmetric 
phase. 

The transition matrix element T(4'' J/^I^ tttt) depends 
at tree level on U only via the mass rn^ of the s me- 
son. We parametrize its temperature dependence by 
(m.(r)2 - m2)/(m0 ' - ml) ^ {{qq){T)/{qq)of\ be- 
cause ms is expected to track roughly the strength of 
the quark condensate (qq). (Other parametrizations are 
also possible, without affecting the main results.) The 
T dependence of the quark condensate has been infered 
from chiral perturbation theory (xPT) calculations to or- 
der T^ with the result that the condensate vanishes 
around 190 (150) MeV for the number of massless flavors 
Nf equal two (three). Several factors limit the applica- 
bility of xPT already below Tc- However, the Tc values 
from this extrapolation are overall consistent with the 
lattice QCD data |jl^. Later we will treat Nf as an ef- 
fective number of degrees of freedom whose magnitude 
controls the value of T^. Above Tc, the scalar and the 
pion mass are taken to be degenerate and T-independent, 
since only a narrow window of temperatures above Tc will 
be probed. 

We flx the vacuum value of the s meson mass to be 1400 
MeV. The particle data group lists a very broad state at 
approximately this mass, the /o(1400). The location of 
the masss of the chiral sigma is controversial in the liter- 
ature. Here we merely note that in the framework of the 
linear sigma model at tree-level a scalar meson with low 
mass is str ong ly disfavored by experimental data on tttt 
scattering [ pTf , |p2| . This still allows a wide range of val- 
ues for the scalar meson mass above 1 GeV. (Note that 
the narrow scalar-isoscalar state just below 1 GeV, the 
/o(975), is probably not the chiral partner of the pions 
but perhaps a KK "molecule".) Threshold theorems are 
not well suited to constrain the s mass much further, be- 
cause the scalar propagator gets effectively replaced by 
a contact term. We determined 5^/47r «0.36 from the 
partial decay width r('I'' 'i'mr). A term im'^Tsim-^Tr) 
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containing the s decay width has been added to the de- 
nominator of the free s meson propagator to take rescat- 
tering corrections into account. Similarly, the s acquires 
a nonzero width from collisions (stt <-> sir) at finite tem- 
perature. 

On the left hand side of Fig. 1 we compare the cross 
section a^'ifiT — > vP'tt) in the vacuum and in the chiral 
phase. Note that the cross section is plotted as a func- 
tion of the out-going pion energy. The in-going pion has 
an energy larger by approximately 600 MeV which is all 
of its typical energy at relevant temperatures. It becomes 
apparent that the thermally averaged cross section is tiny 
if the vacuum values are used. On the right hand side 
the rate constant A = Tded^') + Tco/K*') for transi- 
tions into a J/'^ in the thermal heat bath is displayed. 
rc=170 MeV has been assumed. The total rate shoots 
up at Tc and increases smoothly afterwards. The most 
important contribution to the transition between the two 
charmonium states comes from the collisions with pions. 
For comparison we also show the contribution from the 

decay to the rate and the total rate but with the sigma 
mass at its vacuum value. 
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FIG. 2. Time evolution of *'/* for 3 different initial 
energy densities 1.5/2.5/3 GeV/fm^ and Tc^lTO MeV. The 
equilibrium value at T=Tc is indicated by an arrow. 

Next we would like to study the corrections to the 
^f'/^f ratio, mainly feed-down from x states into J/'^ af- 
ter the collision and non-equlibrium effects. For that we 
consider a simple model which, however, contains the es- 
sential physics. Specifically, we calculate the initial char- 
monium production based on the Glauber model for AA 
collisions at a given impact parameter b. We relate b 
to the experimentally measured transverse energy Et by 
employing the results in Ref. [^. The initial production 
yields of J/\l/, and x in nucleon-nucleon collisions are 
taken as 23, 6.5 and 84.8 respectively. Here the two x 
states are lumped together with average branching ra- 
tio 0.21 into J/'^. The total normalization allows di- 



rect comparison with NA38 and NA50 data referred to 
as Ba{J/'^)/a{DY) in |2|]. A common primordial sup- 
pression of all charmonia depending on parameter L (see 
is assumed as in We calculate the final 

yields by integrating the coupled rate equations for the 
local densities of the charmonium states, e.g. for the 

dN'/dT^~\{^' ^^){N' -N'^q)^\^-,{^')N' . (2) 

Boost invariance is assumed with formation time 
ro=2fm/c. As motivated in Q the rate constants 
Acc for charmonium dissociation are taken as zero be- 
low and infinity above some critical energy density Ccr- 
The ^' equilibrium density iV^^ is related to the J/^ 

density N via N'^^^N{M' /Mf/'^ex:p{-{M' - M)/T). 
Thus we need to model the temperature evolution of 
the created medium. The initial energy densities are 
taken as proportional to the density of participants 
in the transverse plane as in more precisely e = 
dAA{Et/Np)/TodNp/cfr. The choices dsu=2-2 and 
rfpbPb=4.0 lead to maximum local energy densities of 
2.3 and 3.5 GeV/fm'^ for the most central S-l-U and 
Pb-|-Pb collisions. The equation of state which governs 
the isentropic expansion simulates a resonance gas with 
energy density and temperature related by a power law 
e = abT^'^^ and parameters a, 6 taken from j25j. The lat- 
ter choices have been made with an eye on the space-time 
evolution according to RQMD calculations which 
connects our simple model with successful phenomenol- 
ogy of AA collisions. Evolution is stopped at temperature 
140 MeV although the precise value is unimportant due 
to the small transition rates below Tc (cf. Fig. 1). 

The values which determine the energy density 
above which a charmonium state immediately dissolves 
or may not be formed at all are the remaining param- 
eters of the model. They have been adjusted to repro- 
duce the centrality dependence of and J/'^ yields in 
S-FU and Pb-f-Pb collisions (ecr=3.2/l. 5/2.4 GeV/fm^ 
for ^f/^Vx). Fig. 2 illustrates the effect of ecr(*') by 
displaying the time evolution of for different initial 

conditions. At initial energy density just below the crit- 
ical value of 1.5 GeW /iw? the \l/'/\l/ ratio drops sharply 
until the system cools down to the critical temperature 
(here Tc=170 MeV). However, the initial yield is so 
far off equilibrium that the vj/' — > xj/ transitions do not 
equilibrate the yield. What tcr effectively does is to de- 
fine rather crudely a core and a corona region. The larger 
than equilibrium values of ^''/^' in the corona result in 
increasing ^''/^E' values for more peripheral reactions. In 
contrast, starting with larger energy densities above Ccr 
and thus zero initial yield we see that the reactions 
equilibrate the ^' yield. So, /'^ works well as a ther- 
mometer, provided the initial energy densities in colli- 
sions are clearly above 2 GeV/fm'^. 

Fig. 3 contains a direct comparison of the calculated 
and J/* yields for S-l-U and Pb-|-Pb collisions with 
the NA38 and NA50 data. ^' /"^ is displayed for three 
different choices of (150, 160 and 170 MeV). We find 
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excellent agreement with the data for Tc=170 MeV. The 
J/'^ yield is only very mildly affected by the Tc varia- 
tions. The strong thermal ^' suppression factor N^^/N 
works against the vj/ ^ v[/' transitions to be effective as 
a destruction mechanism for J/^'s. On the other side, 
we provide an explanation why scaling from p+A works 
so amazingly well for the J/^ in S+U although some 
of the final feed-down from is gone. As noted above 
the feeds the J/^ yield in the corona of the reaction. 
Thus suppression is partially offset by enhancement 
of direct J/'^'sl 




20 40 60 80 
Et (GeV) 



50 1 00 1 50 
Et (GeV) 



FIG. 3. Calculated J/*]/ yield - normalized to Drell-Yan 
cross section - times branching ratio B into i-i'^fi~ versus pre- 
liminary data from NA38 and NA50 for S on U and Pb on 
Pb (upper panel), 'I''/*!' ratio where each charmonium state is 
weighted with its decay probability into a dimuon pair versus 
data (lower panel). The calculated 'i!' /'i/ values are shown for 
three choices of Tc'. 170 (straight line), 160 (dashed line) and 
150 MeV (dotted line). 

Let us shortly comment on feed-down from x states. 
The feed-down corrections are sizable, even in central 
Pb on Pb collisions. The ^f'/^ ratio changes from 5.64 
% before feed-down to 4.07 % after feed-down. We see 
that both surface emission and feed-down affect ^''/^E' by 
approximately 40 %. However, they contribute with dif- 
ferent sign and cancel almost each other. Furthermore, 
Tc is rather well determined from even if one takes 

the feed-down corrections as a measure of the uncertainty 
involved. To summarize, we have suggested that ^f'/^ 
levels off in ultrarelativistic nucleus-nucleus collisions, be- 
cause 4*' is re-created from J/'?. Transitions between "if' 
and J/'? equilibrate the ratio around Tc due to a drop- 
ping sigma meson mass. A chiral symmetry restoration 
temperature close to 170 MeV explains nicely why ^f'/^* 
approaches a 'universal' value of around 4 %. The future 
dilepton experiments at RHIC and LHC will be impor- 
tant to test the suggested universahty of the final 
ratio. 
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